Effect of Coulomb interactions on the optical properties of doped graphene 
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Recent optical conductivity experiments of doped graphene in the infrared regime reveal a strong 
background in the energy region between the intra and interband transitions difficult to explain 
within conventional pictures. We propose a phenomenological model taking into account the 
marginal Fermi liquid nature of the quasiparticles in graphene near the neutrality point that can 
explain qualitatively the observed features. We also study the electronic Raman signal and suggest 
that it will also be anomalous. 



PACS numbers: 

Introduction. The role of many body corrections to the 
physics of graphene is at this point uncertain. While in 
the first transport experiments electron-electron interac- 
tions seemed not to play a major role [J [1, Q , more recent 
measurements indicate the possible importance of many 
body corrections to the intrinsic properties of graphene 
0, [E H 0, H, H, O- The electronic transport of sin- 
gle layer graphene is one of the most intriguing aspects 
that still lacks a full understanding. Unlike common 
Fermi liquids, graphene shows a minimal conductivity 
at zero frequency in the neutral suspended samples with 
a value of the order of the conductivity quantum, e 2 /h 
that grows linearl y w ith the electron density in the low 
density regime pi. fill. [l2T] . 

The optical properties of graphene (see [Hj] for a very 
accurate reviewVhave been the focus of interest of recent 
experiments H, [1(1 G3 ■ In the visible region they al- 
low a very precise determination of the value of the fine 
structure constant. The structure of graphene is similar 
to that of a two-dimensional two band semiconductor. 
The dynamical conductivity is expected to have a char- 
acteristic shape where intraband (near zero energy) and 
interband (at an energy of twice the chemical potential, 
2/x) transitions can be clearly identified. Recent observa- 
tions in the infrared region of the spectrum [l(|, reveal 
a substantial background of approximately 30 percent of 
the saturation value in the frequency range between in- 
traband and interband transitions. The observed back- 
ground remains approximately constant for several val- 
ues of the chemical potential that range from 0.10 eV up 
to 0.30 eV. Moreover, the threshold at 2/j, shows a very 
wide width independent of gate voltage and therefore of 
carrier density. 

A strong background was also observed in the early 
times of high-T c superconductors and prompted the 
marginal Fermi liquid ideas (lo| . Together with the ob- 
servation of the anomalous infrared optical conductivity 
in [Io| . they observe an increase of the Fermi velocity 
with lowering the frequency at low energies which is sug- 
gested to be due to the renormalization of the Fermi ve- 
locity from the electron-electron interactions as predicted 
in early works [l(| [13] • The measures of optical proper- 
ties as a function of the frequency for varying chemical 



potential 0, [3], angle resolved spectroscopy (ARPES) 
[EH Ell or scanning tunneling microscopy (STM) Q are 
complementary probes to explore the influence of disor- 
der and interactions on the properties of graphene. In 
this work we propose that the mid infrared residual con- 
ductivity found in the experiments at finite chemical po- 
tential can be a footprint of the anomalous behavior of 
the electron lifetime at the Dirac point [l9(. We com- 
pute the optical conductivity with a phenomenological 
self-energy based on the marginal Fermi liquid behav- 
ior discussed in [lj| and see that it helps to understand 
the experiments. We predict a similar anomalous back- 
ground in the electronic contribution to the Raman signal 
which has not been measured yet in graphene. The influ- 
ence of electron-electron interactions on the conductivity 
of undoped graphene was addressed in [2(| HH • 

Main electronic properties of graphene. As it is well 
known by now, graphene is a two dimensional array of 
carbon atoms with an sp 2 coordination. An updated re- 
view of its properties can be found in (22j . The strong sp 2 
sigma bonds form a honeycomb lattice and the low energy 
electronic structure comes from the tt electrons perpen- 
dicular to the plane located at the sites of the honeycomb 
structure. Due to the special topology of the honeycomb 
lattice that has two atoms per unit cell and lacks inver- 
sion symmetry, the Fermi surface of the neutral system 
consists of two inequivalent Fermi points (23| and the 
low energy excitations around them are described in the 
continuum limit by the two dimensional massless Dirac 
Hamiltonian [2~jj . The low energy dispersion relation of 
the two bands around a Fermi point is 

e(k)=±w F |k|, 

from where it can be seen that the density of states van- 
ishes at the Fermi energy and the effective mass of the 
quasiparticles is zero, vf is the Fermi velocity estimated 
to be v F ~ 10 6 m/s ~ c/300. 

The low energy nature of the interacting system was 
studied with renormalization group techniques in [l6| . 
Due to the vanishing of the density of states at the 
Fermi energy the Coulomb interaction is unscreened 
and the neutral system departs from the usual Fermi 
liquid behavior. In particular the inverse lifetime of 
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FIG. 1: (Color online) (a) Measured optical conductivity of graphene for different gate voltages (from right to left V g — 71, 54, 
40, 28, 17 and 10 V). The dashed red line is the optical conductivity at zero chemical potential, (b) Conductivity computed 
with the self-energy (|3]) for parameters a — 0.2, b = 0.001 eV, A = 0.15 eV and for the same gate voltages. The vertical dashed 
lines indicate the 2Ef threshold at each chemical potential. 



the quasiparticles was found in [19j to decrease linearly 
with the frequency instead of the squared Fermi liquid 
behavior. 



around the zero doping and chemical potential. We will 
take this fact into account by considering the following 
phcnomenological electron self-energy: 



Phenomenological model for the electron self-energy. 
The main idea of this work is that the physics of the Dirac 
point influences the infrared properties of the doped sys- 
tem. This assumption has some experimental support. 
In the photoemission analysis done in [25| it is shown 
that interactions deform the spectral function of doped 
graphene not only at the Fermi level but also at the posi- 
tion of the Dirac cone. In the early work [l9| it was shown 
that the inverse lifetime of intrinsic graphene grows lin- 
early with energy at low energies, an infrared behavior 
typical of the marginal Fermi liquid (l5|. The inverse 
lifetime of electrons in the presence of the unscreened 
Coulomb interaction was found to be: 
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This result has been corroborated by recent experiments 
0| that provide an estimation of the coefficient of the 
linear scattering rate of the order of 0.07 in a range of 
energies up to ~ 0.15 eV [34]. 

From a linearly dependent imaginary part of the elec- 
tron self-energy Im E(o>) = a\w\, and computing the real 
part by a Kramers-Kroning transformation we get the 
"marginal Fermi liquid" type of the electron self-energy: 
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where we have also added a constant scattering rate. 
While in the physics of cuprates the electron self-energy 
shows a marginal Fermi liquid behavior for a large range 
of doping levels 15], in the case of graphene the marginal 
behavior is restricted to a very small energy region 
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where A is a low energy cutoff that will be set to A ~ 0.15 
eV. As mentioned before realistic values for the parame- 
ters a, b estimated from the experiment 0] are a = 0.07, 
b = 0.01 eV although they may vary from sample to sam- 
ple. 

A Kramers-Kroning transformation of ([3]) gives for the 
real part 
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Optical properties of the model. From the Kubo for- 
mula we can write the optical conductivity as a function 
of the electron spectral function A s (k,u>) in the so-called 
bubble approximation that neglects the vertex correc- 
tions El: 



Re a(ui) = g s g v 
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J de[72 F {e)-n F (e + uj)]A s (k,Lu)A s ,(k,e + Lj), (5) 

where g StV are the spin and valley degeneracies, Uf(c) = 
[exp(/3(e — fi)) + 1] 1 is the Fermi distribution function, 
and A s (k, to) is the spectral function: 
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FIG. 2: (Color online) (a) Comparison of the optical conductivity computed with the self-energy ([2)| for parameters a = 0.2, 
b = 0.001 eV, A = 0.15 eV, fi — 0.18 eV (red solid line) and same with a constant scattering rate only i. e. a — (blue 
dashed line), (b) Same comparison for the Raman signal. The inset shows the evolution of the Raman signal obtained with 
the self-energy <(3j for decreasing values of the chemical potential from top to bottom. The same color code as in the optical 
conductivity is used. 



The index s = ±1 refers to the upper (lower) part of 
the dispersion relation. The bubble approximation works 
well if the electron self-energy is not strongly k-dependent 
(2?J as happens in our case. 

Finally for systems with a very isotropic Fermi surface 
as it is our case the electronic contribution to the Raman 
signal Im x( w ) can De related to the optical conductivity 
by the simple expression (28j : 

w Re <t(w) ~ Imx(w). (7) 

Results and discussion. We have computed the optical 
conductivity from eqs. and ©. 

Fig. [T]is the main result of this work. It shows a com- 
parison between the dynamical conductivity measured in 
ref. [l(| (a) and that obtained with the self-energy given 
by eq. ([3]) (b) for the following values of the parameters: 
a = 0.2, b = 0.001 eV, A = 0.15 eV. In the experimental 
figure it can be seen a contribution coming from de Drude 
peak at low energies, the interband transitions beyond 2/i 
and the strong background between the Drude peak and 
the interband transitions. It can also be appreciated that 
the width of the threshold of the interband transitions at 
2/i is very wide and almost independent of the electronic 
density. In our approach we obtain the most prominent 
characteristics of the experimental results as can be ap- 
preciated in Fig. DJb). 

In order to clarify the role of the Coulomb interactions 
encoded in parameter a of ([3]) versus a simple constant 
scattering rate b, we have shown in Fig. [2] (a) a com- 
parison of the optical conductivity computed with the 
self-energy @ for parameters a — 0.2, b — 0.001 eV, 
A = 0.15 eV, fi — 0.18 eV (red solid line) and the same 
with only a constant scattering rate (a=0) (blue dashed 
line). For the constant scattering rate the intermediate 
frequency region between Drude and interband transi- 
tions is almost zero and the threshold is sharp. We have 



checked that the same kind of behavior of the conduc- 
tivity is found using a Fermi liquid type of self-en ergy 
which is expected to correspond to doped graphene [2!| • 
When the linear contribution to the scattering rate is 
added the Drude peak becomes wider, there is a consid- 
erable background at intermediate frequencies and the 
threshold of the interband transitions acquires a finite 
width. All these characteristics are expected from a lin- 
ear scattering rate since it induces a general shift of the 
spectral weight from lower to higher frequencies. It is 
worth mentioning that although we have added the cut- 
off A in expression @ to restrict the marginal behavior 
to low frequencies, the shift of the spectral weight is quite 
remarkable. 

The problem of the anomalous residual conductiv- 
ity was studied in [3(J, HH taking into account unitary 
scatterers, charge impurities and optical and acoustic 
phonons. They arrive to the result that to account for 
the observed value of the residual conductivity another 
mechanism effective at high energies must be included. 
They also obtain that optical phonons give a marginal 
contribution similar to Coulomb interactions for frequen- 
cies away from the Dirac point. Some types of disorder 
also induce a linear behavior on the scattering rate at low 
energies [32j . 

What we see from our analysis is that including a linear 
imaginary part of the self-energy even for a very low en- 
ergy range stabilizes the background of the conductivity 
and widens the threshold of the interband transitions in a 
way that is almost independent of the electronic density. 
To get a better fit for the experimental data we have 
chosen a big value of the slope of the linear scattering 
rate "a". We must mention that besides the electronic 
interaction there can be other contributions to the lin- 
ear scatterin g ra te coming from disorder [32| or optical 
phonons [3fj|.l3lT|. 
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Fig. [2jb) shows the same comparison for the electronic 
contribution to the Raman signal in graphcnc yet to be 
measured. Dashed blue curve refers to the case with a 
constant scattering rate. The electronic Raman scatter- 
ing has been calculated very recently in Ref. [33|]. In 
agreement with our work they also find a linear spectral 
density beyond 2/i coming from interband transitions. 
We obtain a Drude contribution at low energy coming 
from the constant scattering rate b = 0.001 eV (dotted 
blue line). As happens in the optical conductivity there 
is no spectral weight in the region between intra and in- 
terband transition when only a constant scattering rate 
is considered. Adding the marginal like behavior however 
(red solid line) there is a noticeable background at inter- 
mediate frequencies and there are no sharp separations 
between the Drude contribution, the background and the 
interband transitions. The inset of Fig. [UR) shows the 
variation of the Raman spectrum with the chemical po- 
tential. The values plotted are the same as in the optical 
conductivity. The three main characteristics i.e. : Drude 
contribution, interband transitions and background are 
clearly seen. For lower chemical potentials it can be seen 
that the Drude peak shifts to higher frequencies when 
increasing fi since the linear contribution dominates the 
physics. This behavior will also be present in the Drude 



contribution of the optical conductivity but such low fre- 
quencies are not accessible in optical conductivity exper- 
iments. 

To summarize, we have proposed a phenomenological 
model that includes the linear behavior of the scattering 
rate at low frequencies and shown that it can explain 
qualitatively recent experiments in optical conductivity 
of graphene The linear scattering rate might 

have its origin on electronic interactions as discussed in 
this paper but optical phonons and disorder can also 
contribute to strengthen the marginal Fermi liquid type 
of behavior. In particular, we have obtained the strong 
background observed between the Drude peak and inter- 
band transitions and the wide width of the threshold of 
the interband transitions, two features that are difficult 
to explain using a Fermi liquid picture. We have pre- 
dicted that an anomalous background will be also seen in 
other optical probes such as electronic Raman scattering. 
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